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We investigate the necessary and sufficient condition for a convex cone of positive semidefinite 
operators to be fixed by a unital quantum operation (j> acting on finite-dimensional quantum states. 
By reducing this problem to the problem of simultaneous diagonalization of the Kraus operators 
associated with <f>, we can completely characterize the kind of quantum states that are fixed by <f>. 
Our work has several applications. It gives a simple proof of the structural characterization of a 
unital quantum operation that acts on finite-dimensional quantum states — a result not explicitly 
mentioned in earlier studies. It also provides a necessary and sufficient condition for what kind 
of measurement statistics is preserved by a unital quantum operation. Finally, our result clarifies 
and extends the work of St0rmer by giving a proof of a reduction theorem on the unassisted and 
entanglement-assisted classical capacities, coherent information, and minimal output Renyi entropy 
of a unital channel acting on finite-dimensional quantum state. 
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I. INTRODUCTION 

A quantum channel can be modeled by a quantum op- 
eration on a separable Hilbert space T-L. Mathemati- 
cally, 4>: B{ri) — > B(H) is a trace-preserving completely 
positive map on the set of all bounded operators of T-L. 
Characterizing quantum operations and studying their 
properties are two important areas of research in quan- 
tum information science. Nevertheless, many apparently 
simple and basic questions regarding a quantum opera- 
tion are difficult to answer. One example is to find all 
quantum states p (in other words, trace one positive self- 
adjoint operators) that are fixed points of a given quan- 
tum operation in the sense that <f>(p) — p. This question 
is still open. Recently, some progress has been made 
on attacking this question. Kribs [_J, who used simple 
functional analysis argument, and Arias et al. 0, who 
applied the generalized Liiders theorem, independently 
discovered a useful necessary and sufficient condition for 
a class of quantum channels known as unital quantum 
operations to fix a quantum state provided that the di- 
mension of the Hilbert space T-L is finite. Studies on the 
generalization and limitations of Arias et al.'s result along 
the line of generalized Liiders theorem have also been re- 
ported. jj_-ti(J 

Restricting the study to unital rather than general 
quantum channels is a sensible tactic to make progress. 
This is because many physical processes in actual ex- 
periments such as depolarization and dephasing can be 
modeled by unital channels so that it is worthwhile to 
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study these channels. Besides, the mathematical tools 
and results to deal with unital quantum operations are 
reasonably well-developed. Therefore, it is not surpris- 
ing that several quantum information science problems 
involving unital quantum operations have been solved. 
For example, by means of finite-dimensional C*-algebra, 
St0rmer found that the evaluation of the (unassisted) 
classical capacity Ci j00 of a unital quantum channel act- 
ing on finite-dimensional quantum states can be reduced 
to calculating the same channel capacity for the case in 
which the image of a non-scalar projection is never a 
projection. [ll| In another study, Blume-Kohout et al. 
used unital quantum channel as an auxiliary tool and re- 
sults from matrix algebra to characterize the geometric 
structure of noiseless subsystems, decoherence-free sub- 
spaces, pointer bases and quantum error-correcting codes 
acting on finite-dimensional quantum states. The geo- 
metric structure they have identified is an isometry to 
fixed points of certain unital quantum operations. [l2.[l3| 
They also extended their results by showing several inter- 
esting properties of the set of fixed points of an arbitrary 
finite-dimensional quantum channel. [l__[ Using similar 
techniques, Rosmanis studied the properties of the fixed 
space of a positive (but not necessarily completely pos- 
itive) trace-preserving map. [___| Recently, Mendl and 
Wolf discovered several equivalent definitions for a unital 
quantum operation and used them to study the relation 
between unital channels and quantum error corrections. 
Their results shine new light on the asymptotic Birkhoff 
conjecture and the separation of the set of mixtures of 
unitary channels. (To| Lately, Zhang and Wu showed an 
easily checkable necessary and sufficient condition for a 
finite-dimensional quantum state whose von Neumann 
entropy is preserved by a given unital channel. Interest- 
ingly, their work is closely related to finding fixed points 
of the composition of a unital quan tum operation <f> and 
its adjoint quantum operation. [16j 
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In turns out that the works of Kribs 0, St0rmer [n| 
and Blume-Kohout et al. [HI, [l3[ mentioned in the above 
paragraph are closely related to the following variation 
of the fixed point problem in which we called the fixed 
convex cone of positive semidcfinitc operator problem. 
Let B(H) be the set of all bounded operators of a finite- 
dimensional Hilbert space H. (All Hilbert spaces con- 
sidered in this paper are of finite dimensions.) And let 
4>: B{Ti) — I B(H) be a unital quantum operation. That is 
to say, 4>(p) can always be expressed as a finite sum in the 
form AipA\ known as the operator-sum representa- 
tion with Y,i A** = 7 « = E 2 MA so that ^( 7 w) = J «> 
where 1% denotes the identity operator on H. Surely, A^s 
(which are called the Kraus operators associated with 0) 
as well as elements in B{%) can be regarded as complex- 
valued matrices of dimension dim % < oo with respect 
to an orthonormal basis of H. By the fixed convex cone 
of positive semidefinite operator problem, we refer to the 
problem of determining if (f> fixes a convex cone C + formed 
by the set of all positive semidefinite operators in B(S) 
for a (proper) Hilbert subspace S of T-L in the sense that 
</>[C + ] C C + . Since </) is a trace-preserving completely pos- 
itive map, our fixed convex cone problem is equivalent to 
the problem of determining the existence of a subspace S 
of % such that cj> sends density matrices in S to density 
matrices in S — a problem of quantum information sci- 
ence interest. (Unlike the approaches in Refs. (ll [Tl| - [l3| . 
we are not interested in the action of <j> on a e B(TL) that 
is not a density matrix because a docs not describe a 
physical quantum state and hence plays no role in quan- 
tum information science.) By abusing language, we call 
the Hilbert subspace S an invariant subspace of <fr; alter- 
natively, we say that </> fixes the subspace S. If such a 
proper subspace S exists, we would like to explicitly find 
it out. 

Interestingly, we are going to show in Sec. |H] that if 
S is an invariant subspace of a unital quantum opera- 
tion <fi, then so is its orthogonal subspace S . In other 
words, a unital quantum operation <f> induces a direct sum 
decomposition of the finite-dimensional Hilbert space T-L 
into irreducible invariant subspaces (IrlS's) of <fi. Here an 
IrlS means that it does not contain any proper invariant 
subspace. More importantly, we further prove in Sec. |H] 
that 4> fixes the subspace S (and hence also S ) if and 
only if S and S are simultaneous invariant subspaces 
of all Kraus operators Aj's associated with <j>. That is, 
Ai[S] C S and Ai[S ] c S for all i. Hence, <j> has an 
interesting structure in the sense that it induces a simul- 
taneous block diagonalization for all its Kraus operators 
Ai's such that each diagonal block acts on a different 
irreducible invariant subspace of <f>. In this regard, our 
notion of convex cone fixation, which concentrates only 
on the action of <p on density matrices, turns out to be 
strong enough to force <p[B(S)] C B(S). 

An interesting consequence of this finding is a simple 
proof of the structural characterization theorem of uni- 
tal quantum operation acting on finite-dimensional den- 
sity matrices. As far as we know, we are the first group 



who explicitly state and prove this theorem although this 
theorem can be deduced from the works of Kribs [l| and 
Blume-Kohout et al. [l3j . Furthermore, both prior works 
did not investigate the quantum information processing 
consequences of the structural characterization theorem 
of finite-dimensional unital quantum operations. In fact, 
Kribs showed in Lemma 2.2 of Ref. [1| that fixed points 
for an irreducible unital quantum operation acting on 
B(H) with dimH < oo must be scalars. This is a spe- 
cial case of Theorem [2] to be reported in Sec. HU How- 
ever, Kribs did not mention the decomposition of unital 
quantum operation into direct sum of irreducible ones 
and properties of such a decomposition. Whereas in 
Lemmas 5.4, 5.5 and Theorem 5 of Ref. [T3|, Blume- 
Kohout et al. used the structure theorem of matrix al- 
gebra to prove the existence of a direct sum decomposi- 
tion structure for a general finite-dimensional quantum 
operation. They also discussed its implications on the 
structure of the associated Kraus operators. Our struc- 
tural characterization theorem can be readily deduced 
from this work by sharpening their results for the case of 
a unital operation. Here we use an alternative approach, 
which uses rather elementary techniques in mathematical 
analysis and graph theory, to obtain the structural char- 
acterization theorem for finite-dimensional unital quan- 
tum operations. In this way, we avoid going through 
the more technical proofs on a more general situation in 
Refs. [J [l3| and adapting them to our particular case of 
interest. 

The structure theorem reported in Sec. |H] has a few 
quantum information science applications. In Sec. lIIIl we 
first use it to obtain a simple and intuitive proof of a theo- 
rem concerning the calculation of classical capacity Ci j00 
of a unital channel acting on finite-dimensional quantum 



111. Oi 



states originally obtained by St0rmer in Ref. 
proof can be used to extend St0rmer's result to the calcu- 
lation of the entanglement-assisted classical capacity C e , 
the coherent information J, and the minimal output a- 
Renyi entropy S m i n , a of the same channel. More impor- 
tantly, we completely characterize the kind of quantum 
states that are fixed by a unital quantum operation in 
Sec. lIIIl And we provide a necessary and sufficient condi- 
tion for the measurement statistics of a positive operator- 
valued measure (POVM) measurement to be preserved 
when a general finite-dimensional quantum state passes 
through a unital quantum channel. Finally, we briefly 



discuss the implications of our results in Sec. [TV] 



II. STRUCTURAL CHARACTERIZATION OF 
UNITAL QUANTUM OPERATIONS ACTING ON 
FINITE-DIMENSIONAL DENSITY MATRICES 

We use the following theorem, which was indepen- 
dently proven by Kribs as Theorem 2.1 in Ref. [l[ and by 
Arias et al. as Theorem 3.5(a) in Ref. 0], as our starting 
point to study the existence of invariant subspace of a 
unital quantum operation. (See also the variation of this 
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theorem stated as Lemma 5.2 in Ref. [HJ].) 

Theorem 1. Let cft(-) = J^. Ai ■ A\ be a unital quantum 
operation on the set of all bounded operators of a finite- 
dimensional Hilbert space %. Then, <ft fixes a G B(H.) 
(that is, <fi(o~) = a) if and only if Ai<j = o~Ai for all i. 

Remark 1. In particular, Theorem [7] relates fixing a 
quantum state p (a property independent of the choice of 
the associated Kraus operators Ai 's) to the commutativ- 
ity of p with the set of Kraus operators Ai 's used any 
operator-sum representation of the unital quantum oper- 
ator (ft. This is possible partly because of a unitary degree 
of freedom in the operator-sum representation of (ft. More 
precisely, (ft(-) = £\ Ai ■ A\ = £\ &i ' B l and or % */ 
Bj = J2i u ijAi for all j where the (i,j)-th ele- 

ment of a unitary matrix. Thus, the commutativity 
of p with all Ai 's implies the commutativity of p with all 
Bi's. 

A special case of Theorem Q] is that <f> fixes a (normal- 
ized) pure state \x) if and only if \x) is an eigenvector of 
Ai for all i. (The "if part" can be deduced by the fact that 
Ai\x) = \i\x) implies (ft{\x){x\) = J2i W\ 2 \ x )( x \ = \ x )( x \ 
for (ft is trace-preserving. The "only if part" follows from 

Ai\x) = Ai\x){x\x) = \x){x\Ai\x) = (x\Ai\x) \x). (1) 

Note also that by the same argument in Remark [TJ |a:) 
is an eigenvector for any Kraus operators used in the 
operator-sum representation of (ft although the corre- 
sponding eigenvalues are operator-sum representation de- 
pendent.) In other words, (ft fixes an one-dimensional sub- 
space S of H if and only if S is a simultaneous eigenspace 
of all its Kraus operators A^s. 

We now apply Theorem [T] to study the necessary and 
sufficient condition for the existence of invariant subspace 
of 0. 

Theorem 2. Let (ft be a unital quantum operation on 
the set of all bounded operators of a finite- dimensional 
Hilbert space %. Then, (ft fixes a Hilbert subspace S of 
TL if and only if (ft(Ps) = Ps where Ps is the projection 
operator onto S. 

Proof. {<==)'■ Assume there exists an Hermitian operator 
as whose support is in S. Suppose further that < 
fT 5 < Ps and (ft(as) & S. Then, a s > = Ps - os > 0. 
More importantly, 

Ps = 4>(Ps) = 0(o-s ) + <j>{o- s -)- (2) 

Since <p(os) & S, there exists |x) G S such that 
(x\(f>(as)\x) ^ where S 1 - denotes the orthogonal com- 
plement of S. And by the positivity of <p(o~s), we have 
(x\(f)(as)\x) > 0. But then 

(x\(j)(asi)\x) = (x\P s \x) - (x\<p(os)\x) 

= -(x\4>(vs)\x)<0. (3) 



This is impossible for is a quantum operation. There- 
fore, we conclude that 4>{o~s) G S for all non- negative 
operators on S. In other words, (ft fixes S. 

(=*-): 4> fixes S implies 4>(Ps) £ S. Assume the con- 
trary is true so that <f>(Ps) ^ Ps- Let P s ± = Pu — Ps be 
the projection operator onto the orthogonal complement 
of S. Then, 

(ft{P s ) + - <KPh) = Pu- (4) 

Note that the last equality in the above equation fol- 
lows from the fact that (ft is unital. Since Ti is finite- 
dimensional and (ft fixes S, we can express (ft(Ps) as the fi- 
nite sum X^™ 5 a j\ v j)(yj \ w hh non-negative aj's, where 
{!%)} i s an orthonormal basis of S. Since (ft(Ps) 7^ Ps, 
we cannot have aj = 1 for j = 1,2,..., dim S. Neverthe- 
less, since (ft is trace-preserving, Oj's still have to satisfy 

the constraint Y^j=l S a j = dim 5. Thus, by relabeling 
the index if necessary, we may assume that a\ > 1 and 
a 2 < 1. So, from Eq. g]), 

1 = ( yi \<f>(P s )\yi) + (yi\<KP s ±)\vi) 
= a 1 + (y 1 \4>(P s ±)\y 1 ), (5) 

That is to say, {yi\<ft{Ps ± )\yi) = 1 — oi < 0, which con- 
tradicts the assumption that is a quantum operation. 
Therefore, we conclude that <ft(Ps) = Ps- D 

Remark 2. Note that the unital condition is needed only 
in the proof of the "only if part" of the Theorem. Note 
further that Kribs proved the special case of this theorem 
when (ft does not fix any proper subspace of H Jj/. 

Corollary 1. (ft fixes S if and only if (ft fixes S . 

Proof. By Theorem [5] and Eq. (QJ, (ft fixes S if and only if 
4>{P S ± )=Ph- 0{Ps) =Pu-Ps = P S ^ , which in turn 
is true if and only if (ft fixes S . □ 

Corollary 2. (ft fixes the subspace S if and only if 

(s x \A i \s)=Q = (s\A i \s x ) (6) 

for all \s) G S and |s^} G S 1 - and for all Kraus operators 
Ai 's of (j). 

Proof. From Theorem [2j (ft fixes S implies (ft(Ps) = Ps- 
Theorem [T] further implies AiPs = PsAi for all i. 
Multiplying (s | on the left and \s) on the right gives 
(s |A,-|s) = 0. Similarly, multiplying (s\ on the left and 
Is- 1 ) on the right gives (s|^4i Is^) = 0. 

To prove the converse, one only needs to observe that 
Eq. © implies Ai[S] C S for all i. Thus, <j) fixes S. □ 

Remark 3. Surprisingly, the notion of convex cone fix- 
ation by a unital quantum operation (f> is much stronger 
than what we have originally written down. Recall from 
Sec. that (ft fixes a Hilbert subspace S simply means 
<ft[Cg] G Cg where C$ denotes the convex cone formed 
by the set of all positive semidefinite operators in B(S). 
Yet, the above two Corollaries say that subspace fixing 
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actually demands something more, namely, <p[B(S)] C 
B(S), cj)[B(S^)\ C B(S ± ). Note further that we cannot 
directly deduce Corollary [H from Corollary [7] for the lat- 
ter says nothing on 4>(o~) for a general operator a whose 
support is in S. 

This motivates us to formulate our main theorem. 

Theorem 3 (Structural theorem for unital quan- 
tum operations on finite-dimensional density matrices). 

(a) Every finite- dimensional unital quantum operation 
4>: B(H) — > B(H) induces a direct sum decomposition of 
H = Sj where each Sj is an IrlS of <j). Furthermore, 
every Kraus operator Ai of (f> can also be decomposed as 

A i 3 where A i 3 € B(Sj) for all i,j. Hence, the quan- 
tum operation 4>\ B ^ S j is also unital whenever S is an in- 
variant subspace of (p. In other words, <f> can be expressed 
as a direct sum^ - 4>\ B ^s-) of unital quantum operations. 
In matrix language, <j) has a proper invariant subspace if 
and only if each of its Kraus operators can be simulta- 
neously block diagonalized by unitary conjugation into at 
least two diagonal blocks. 

(b) In addition, S is an IrlS of <f> if and only if all 
fixed positive self-adjoint operators of <f> in B(S) are in 
the form aP$ for some a > 0, where Ps denotes the 
projection operator onto S. 

Proof. From Corollaries Q] and [2j every Kraus operator 
Ai of a unital quantum operator </>: B(H) — > B(H) with 
dimH < oo admits a direct sum decomposition Af(BAf 
where Af (Af ) is a linear operator in the invariant sub- 
space S (S^) of <fi. Since J2i A f Af^ = Is, the quantum 
operation 4>\ B ( S ) (•) = J2i Af-Af^ is unital. By the same 
token, 0|g( S ±) 

(') = A f ■ A f is unital. By recur- 
sively applying Corollaries [1] and [2] to the unital quantum 
operations ^Igrj) (•) and <^>|g(5-M at most dim (H) — 1 
times, part (a) is proven. 

To prove part (b), suppose < a G B(S) is fixed by <j> 
and yet a ^ aPs for all a > 0. Since T-L and hence S are 
finite-dimensional, we may write a — Ej™ 5 
for some orthonormal basis vectors ) 's of S and all 
fej's are non-negative. By relabeling the indices if nec- 
essary, may we assume that b\ < bj for all j; and the 
requirement that a ^ aPs implies not all bj's are equal. 
Using part (a), (f>\ s is also unital so that tfi fixes Ps and 
hence also the operator a' given by 

dim S 

a' = a- b.Ps = ]T (b s - h) \Vj)(Vj\ > 0. (7) 

From Theorem [TJ Af a' = a'Af for all i. By multi- 
plying {yi\ on the left and \yj) on the right, we get 
(bj — b{)(yi\Af \yj) = for all i,j. Similarly, mul- 
tiplying (yj\ on the left and \yx) on the right gives 
(bj -h){yj\Af lift) = for all That is, (yj\Af\ yi ) = 
(yi\Af\yj) = for all i whenever b\ < bj. Following the 



same logic, we conclude that 

(y j \Af\y k ) = (y k \Af\y j )=0 (8) 

for all i whenever b\ = bk < bj . Consequently, each of the 
Af 's can be simultaneously block diagonalized by unitary 
conjugation to at least two diagonal blocks — one block 
corresponds to those indices fc's with bk — b\ and the 
other block corresponds to those fc's with b^> b\. From 
Corollary S is not an MS of T-L, which is absurd. 

Finally, to show the converse of part (b), suppose S' is 
a proper subspace of S that is fixed by <f). Then, it is clear 
from part (a) that (j) fixes the projection operator P$' ■ 
Hence, not all positive self-adjoint operators in B(S) fixed 
by <p are in the form aPs. This proves the theorem. □ 

Remark 4. As we have mentioned in Secjft Theorem^ 
can also be deduced from the work of Kribs fj/ and Blume- 
Kohout et al. fp\] . 

Remark 5. The unitary degree of freedom in operator- 
sum representation mentioned in Remark[]]is the reason 
why Eq. ([6]) in Corollary [H as well as the simultaneous 
block diagonalization of Ai's in part (a) of Theorem^ 
hold for any Kraus operators associated with <fi. 

Remark 6. Actually, positivity of the self-adjoint oper- 
ator a is not essential in the proof of part (b ) of Theo- 
rem [3J In fact, it can be slightly strengthened as S is an 
IrlS of (j) if and only if all fixed self-adjoint operators of 
(j) in B(H) are in the form aP$ for some a G M. The 
proof is left to interested readers. We do not state this 
slightly more general form in the theorem because we are 
more interested in the action of (j) to a physical quantum 
state. 

Remark 7. In general, a matrix can be diagonalized into 
irreducible blocks in more than one orthonormal basis due 
to degeneracy of its eigenspace. This is also the reason 
why the IrlS decomposition of a finite- dimensional unital 
quantum operation discussed in the above theorem need 
not be unique. 

The following example illustrates the power of Theo- 
rem [3] in determining the structures of some well-known 
unital quantum channels. 

Example 1. Consider the depolarization qudit channel 
<f> over a finite- dimensional Hilbert space %. That is, 
<f> sends a density matrix p to (1 — p)p + pin / dim % 
with < p < 1. More generally, (f>(o~) = (1 — p)o~ + 
pI-H Tr((j)/dim7^ for all a G B(7i). Let Ps be the pro- 
jector on a Hilbert subspace SofH. Then, cf> fixes Pg if 
and only if pPs = pin dim S/ dim % . As p ^ 0, <f> fixes 
Ps if and only if S = H. Since <fi is unital, Theorem^ 
implies that % is the only IrlS of <fi. In this regard, all 
finite- dimensional depolarizing qudit channels cannot be 
further decomposed as a direct sum of two unital quantum 
operations. 
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Although IrlS decomposition induced by <j> may not be 
unique, different MS decompositions share the same di- 
mensional structure. This is analogous to the dimensions 
of IrlS for complex- valued matrices. 

Theorem 4. Suppose % = 0™ 1 S, = ©? =1 <Sj &e 
<wo different IrlS decompositions induced by q> reported 
in part (a) of Theorem^ indexed in such a way that 
dim Si > dim5j' whenever i > i' and dimSj > dimSj, 
whenever j > j' . Then m = n and dimS, = dimS^ for 
all i. 

Proof. Since there are two distinct IrlS decompositions 
for (p, we can always find an IrlS in the first decompo- 
sition, say, <Si that is distinct from all the IrlS's in the 
second decomposition. As S\ C 0jSj, we can find an 
IrlS in the second decomposition, say S[ which is not 
contained in S^~ the orthogonal subspace of Si. 

We claim that dim Si — dim S[ . Suppose the contrary, 
may we assume without lost of generality that dim Si > 
dim S[ . Let us write Si — Ti ®Ui and S[ = T{®U[ where 
Ti = Si n S'i and T{ =S[n S^. Since dim Si > dimS( 
and S[ <£_ Si , we conclude that 71 is a proper subspace 
of Si. 

Because Si and S[ are IrlS's of <f>, by part (b) of Theo- 
rem [31 cj) fixes the projectors Ps l: Ps[ and hence also the 
positive operator a = Ps x + 0-5Ps[ S B(Si + S[). 

We are going to show that the only eigenvectors of 
a with eigenvalue 1 are vectors in the subspace Ti- A 
vector in Si +S[ can be uniquely written as = a\x) + 
b\x ) where \x) and \x ± ) are normalized vectors in Si and 
(Si + S[) CiSi, respectively. Then \ip) is an eigenvector 
of a with eigenvalue 1 provided that 

(P Sl + 0.5P 5; ) (a\x) + blar 1 )) = a\x) + b^). (9) 

Multiplying (x\ and (x | to Eq. (0) gives the system of 
equations 

a(x\P sl \x) = -b(x\P s >\x ± ), (10a) 
Q.Saix^Ps^x) =b(l-0.5{x ± \P s{ \x ± }) . (10b) 

This system of equations has non-trivial solution if and 
only if 

(a;|P si |a:)(l-0.5(a:- L |Psi|a: x )) = -0.5 \{x x \P s[ \x}\ 2 . 

(11) 

As Pg> is a projector, (y\Ps' \y) G [0, 1] for all normalized 
vector \y). Therefore, the L.H.S. and R.H.S. of Eq. (JTTJ) 
are non-negative and non-positive, respectively. Conse- 
quently, both the R.H.S. and R.H.S. of Eq. (JTTJl equal 0. 
Thus, (x^Ps^x) = 0. Since 1 - 0.5(x ± \P S ' i \x ± ) > 0, 
Eq. (llObp implies 6 = 0. So, = a\x) G 7i is the so- 
lution of Eq. In other words, we have succeeded in 
showing that all eigenvectors of a with eigenvalue 1 are 
the vectors in 7i. 

Therefore, a can be rewritten as Pj\ © a 1 where a' is a 
positive operator whose eigenspectrum does not contain 



the eigenvalue 1 (and hence its support is not in 71). 
Using the argument in the proof of part (b) of TheoremUJ 
we conclude that cj> fixes Pj\ . In other words, 71 is an 
proper invariant subspace of Si. This contradicts the 
fact that Si is irreducible. Therefore, we conclude that 
dim Si = dimS^. 

To summarize, so far we have deduced that if Si <£_ 
S'j (and hence Sj (£_ S^), then dim Si = dimSj. We 
now construct a finite bipartite graph whose vertices 
are the IrlS's Si's and Sj's — each side corresponds to 
a different IrlS decomposition of TL. An edge is drawn 
between vertices Si and Sj each picked from one side 

provided that either (i) Si (jL S'j or (ii) Si = S'-. Our 
construction of © guarantees that: 

• Each IrlS must connect to at least one IrlS on the 
other side of the graph. 

• The dimension of each Hilbert space associated 
with the IrlS's in the same connected component 
of © agrees. 

• The direct sum of IrlS's drawn from each of the 
two sides in a given connect component of © must 
agree. 

Since T-L is finite-dimensional, we further conclude that 
in each connected component, the number of IrlS's in 
the two sides of © are the same. Therefore, in each 
connected component, we can construct a bijection from 
IrlS's in one side to the IrlS's in the other side of the 
graph ©. This bijection maps each IrlS in one of the IrlS 
decomposition to an IrlS of equal dimension in the other 
decomposition. Hence the theorem is proved. □ 

Remark 8. Actually, the construction of the bipartite 
graph and the determination of the connected compo- 
nents of this graph in the proof of Theorem [^] are quite 
efficient. To see this, first by Gram-Schmidt orthogo- 
nalization procedure, we may represent each subspace Si 
of W by its orthogonal basis. Clearly, the computational 
cost needed to specify all the orthogonal bases of Si 's and 
S'j 's, measured in terms of the number of real number 
arithmetical operations, scales like 0((dim%) 2 ). Given 
a normalized vector \v) and a subspace S specified by a 
set of orthonormal basis {\bk}}, we may determine if 
\v) is in S or not by checking if J2i \ ( v \°i)\ equal to 
1 or not within a certain error tolerant level e caused 
by finite-precision arithmetic. Similarly, we know if \v) 
is in S^ or not by checking if J^i \ { v \bi)\ is equal to 
or not within a certain error tolerant level e. Both 
procedures require 0(dimS dim"H) real number arith- 
metical operations. Using these procedures as subrou- 
tines, we can determine if Si (£ S'j' or Si = S'j in 
0(dim Si dimSj dimH) time by checking if each ba- 
sis vector of S is in S'f or S'j. Hence, the bipartite 
graph © described in the proof of Theorem^ can be con- 
structed in^ li j 0(dimS'i dimSj dimH) = 0((dimH) 3 ) 
time. Furthermore, connected components of a graph can 
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be found by the well-known depth first search algorithm 
in graph theory, which needs a computational time lin- 
ear in the number of edges of the graph f/jj /. So, in 
our case, connected components of 25 can be determined 
in at most 0((dim%) 2 ) time. To summarize, the algo- 
rithm concerning the graph © in the proof of Theorem ^ 
takes 0((dim7{) 3 ) time with the graph construction be- 
ing the most time-consuming step. Thus, the algorithm 
of identifying isomorphic IrlS decompositions reported in 
the proof of Theorem ^ is quite efficient and can be used 
in actual situations especially when dimH is large. 



III. SIMPLE APPLICATIONS OF THE 
STRUCTURAL THEOREM IN QUANTUM 
INFORMATION PROCESSING 

We report three simple applications of Theorem [3] that 
are of quantum information science interest. The first 
one concerns a simplified proof and an extension of the 
result by St0rmer in Ref. fll| on the channel capacity 
of a unital channel acting on finite-dimensional quantum 
states. The second one concerns the type of quantum 
states that are invariant after passing through a unital 
quantum channel. The third one concerns the preser- 
vation of POVM measurement statistics under a unital 
quantum channel. Our work here strengthens the find- 
ings of Arias et al. in Ref. 0. 

In this Section, <fi always denotes a unital quantum op- 
eration on B(tH) with Kraus operators A^s and dimH < 



Theorem 5. Let % = Sj be a direct sum decompo- 
sition of the Hilbert space T-L into IrlS's induced by <fi. 
Then 



Proof. By Theorem [3l we know that the can be ex- 
pressed as a direct sum of finite-dimensional quantum 
operations. Then, all the above formulae concerning the 
reduction of various information-theoretic quantities of 
the quantum operation <f> follow from Proposition 1 in 
Ref. [19| by Fukuda and Wolf, which used basic proper- 
ties such as von Neumann entropy and mutual informa- 
tion are concave functions to show that each of the above 
capacities can be achieved by a density matrix respecting 
the direct sum structure of the quantum operation <f>. □ 



Remark 9. Eq. (|i~2dl) was first proven by St0rmer in 
Theorem 3 of Ref. [11] . The rather complicated condi- 
tions in that Theorem is nothing but the structural de- 
composition of a unital 4> according to our Theorem\3(b) . 

Corollary 3. The set of all quantum states that are fixed 
by 4> are the classical mixtures of the completely mixed 
states (dirriiSj) -1 P$j where each Sj is an IrlS of (p. 

Outline proof. We can always write a density matrix 
p E B(H) in the form Ylj^i ^j\Vj){yj\ f° r some or- 
thonormal basis vectors \yi)'s of H with b\ < bj < 
for all j. Suppose 4>{p) — p. Following the same argu- 
ment used in the proof of part (b) of Theorem[3l we know 
that p — b\Pu is a positive operator fixed by (f>. Besides, 
T = span{|yj): bj = 61} and T 1 = span{|z/j): bj ^ 61} 
are two mutually orthogonal invariant subspaces of 4>. As 
H is finite-dimensional, by recursively applying this ar- 
gument to p — biPy, and (f>\j-, a finite number of times 
until p — b\Pu = 0, we conclude that p is a finite sum in 
the form Ylj c o^Tj wnere c j > an d TJ 's are mutually 
orthogonal invariant subspaces of H. Since each Pf. is 
(possibly) a classical mixture of completely mixed states 
in the form (dimiSfc) 1 P$ k 's, this corollary is proved. □ 



the minimal output a-Renyi entropy (a > 1) is 
given by 



S m in,a(4>) = min 
3 



,a (^Ib(S 3 )) 



the coherent information is given by 
J 



J{4>) = max 
3 



16(5/ 



(12a) 



(12b) 



• the entanglement- assisted classical capacity is given 
by 



C e {cj>) =log^2 Cc(0lB ^> ) , 



(12c) 



the (unassisted) classical capacity is given by 



C lt00 {4>)= log 2 



Cl,oo(0| e(S )) 



(12d) 



The above corollary means that completely mixed 
states of IrlS's are the basic building blocks of quantum 
states fixed by <p. In this sense, the problem of find- 
ing and characterizing invariant subspaces of 4> we are 
studying here is more than a variation of the fixed state 
problem of <f>. It is, in fact, a generalization of the fixed 
state problem. 

Remark 10. Actually, the self-adjointness of p is essen- 
tial in the proof of Corollary [21 while the positivity of p 
is not. So, the Corollary can be slightly extended by say- 
ing that all self-adjoint operators that are fixed by 4> must 
be in the form aiPg i with dj G M. Nevertheless, we 
stress that Corollary does not cover the case of fixing 
a non- self- adjoint operator — a situation of no physical 
meaning. 

Remark 11. Using Theorem^ Corollary^ and the fact 
that a matrix admits two non- orthogonal eigenvectors if 
and only if it has a degenerate eigenspace, it is easy to 
see that the following statements are equivalent: 

• <j) fixes two non- orthogonal pure states \x\), \x2). 
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• <f> admits two distinct IrlS decompositions such that 
the first decomposition contains the IrlS generated 
by \x\) and the second decomposition contains the 
IrlS generated by \x2)- 

• \x%), \x%) are degenerate eigenvectors of each of the 
Kraus operators A4 's of cf> and that AiP$ = PsAi 
for all i where Ps denotes the projector onto the 
space spanned by \xi) and \x2). 

• <j) fixes all pure states in the span of \x\), \x2)- 

To address the question of preservation of measure- 
ment statistics, we begin with the following lemma. 

Lemma 1. Let IT be a projector on the Hilbert space H. 
Then 



(13) 



for all p G B(T~l) if and only if the image ofH under II, 
that is, IL[H], is an invariant subspace of (p. 

Proof. If S = H[H] is an invariant subspace of then 
part (a) of Theorem [3J implies that each Kraus operator 
of (j> can be written as Ai = Af © Af using the notation 
of that Theorem. Regarding II and A^s as block matrices 
in a basis compatible with the H = S © direct sum 
structure, 



n<Mp)nt 

=E 

i 

=E 

i 

=E 



Is 




\Af 




\Af 


t 


Is 


s x_ 






P 


. Af\ 







Af 



Af 



AfTlprtAf 



E[ Af 

■4 (npnt) 



Af 





npitf 




\Af 




o 5 ^_ 




. Af\ 



for all p G B(n). 

To prove the converse, we observe that S = H[H] is 
a Hilbert subspace of %. For any density matrix p £ 
B{S), Eq. (HU) becomes II0(»nt = <j>{p). Therefore, 
the density matrix <fr(p) G B(S). So S is an invariant 
subspace of <p by definition. □ 



Theorem 6. (a) Let tp: B(H) — > B(H) be the quan- 
tum operation p H> ^ k rifc/jllj, where lik 's are projectors 
obeying ^2 k Tlk = I-h o-nd 11^11^' = whenever k ^ k' . 
Then 



ipocj){p) = (f)oip{p) 



(15) 



for all p G B{H.) if and only if every IT^ [H] is an invariant 
subspace of (f>. More importantly, 



for all density matrices p in BlTi) and for all k if and 
only ifTlk[H] is an invariant subspace of (f>. 
(b) More generally, consider a POVM measurement on 
a quantum state in the Hilbert space T-L with POVM ele- 
ments {Ek}. (Note that the number of POVM elements 
need not be finite here.) Then 

Tr (E k p) = Tr [E k cj>(p)} (17) 

for all density matrices p in B(TL) if and only if E^ = 
Y]j dj Psj with a,j > where Pg j is the projector on 
an IrlS Sj of <f>. 

Proof. To prove part (a), note that Eq. (fT5j) is a direct 
consequence of Lemma [T] and linearity of <j). 

Suppose is an invariant subspace of <f>. By tak- 

ing trace in both sides of Eq. (ITU)) and by using the 
fact that <j) is trace-preserving, we have Tr [Tlk4>(p)] = 



Tr 



= Tr| 



(UkpUlj =Tr(ILp). 

To prove the converse in part (a), suppose Eq. (fl6|) 
holds. We set p to an arbitrary but fixed density ma- 
trix in B{n k [H}). Then, 1 = Tr {U k p) = Tr [n fe 0(p)n^ . 
Since <p is trace-preserving and positive, IL is a projector 
and dim"H is finite, we conclude that cf>(p) G B(Hk[H]). 
Hence, nfe[%] is an invariant subspace of 4>. 

We now move on to prove part (b). The sufficiency 
condition in part (b) is a direct consequence of part (a) 
which says that Tr (P Sj p) = Tr [P Sj( j)(p)] for all IrlS Sj. 

To prove the condition is necessary in part (b), we 
use the fact that each Ek is a positive operator in BCH) 
and hence self-adjoint. As % is finite-dimensional, Ek 
can be written as the finite sum y~] ( bjP-r, where 7^'s 
are mutually orthogonal subspaces of H, be > and 
bg < be/ whenever I > £'. (Surely, 6/s and TVs de- 
pend on k. But we do not explicitly emphasize this 
dependence to avoid clumsy notations.) Consequently, 
Tr (Ekp) < bi for all density matrices p G B(7i) with 
equality holds if and only if p G B{T\). Let pi G B(Ti). 
Eq. (HZJ implies Tr [E k (f>{pi)] = b x . Hence, (f>{pi) G B(Ti) 
and 7i is an invariant subspace of ip. By Theorem [31 
( 14 ) 0[B(71 X )] C BiT^). More importantly, Tr [E k p'] < b 2 
for all density matrices p' G S(7i ) with equality holds 
if and only if p' G Biffe)- So, by inductively applying the 
previous argument, we conclude that all 7^'s are invariant 
subspaces of 4>. This proves the converse in part (b). □ 



Tr (U k p) = Tr p fc 0(p)] 



(16) 



Remark 12. The above Theorem says that the order 
of passing a quantum state through the unital channel <p 
and performing a projective measurement ip on the state 
is not important provided that dimH < 00. Moreover, let 
us consider the following two machines — the first one 
performs a POVM measurement on an input quantum 
state and the second one performs the same POVM mea- 
surement after passing the input quantum state through a 
unital quantum channel. The above Theorem completely 
characterizes the kind of POVM measurements in the 
above two machines so that they have the same measure- 
ment statistics given any input quantum state. One im- 
plication of our findings is that if we do not have control 
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on the kind of input quantum states, then measurement 
statistics may be changed by a unital quantum channel <p 
if the measurement is finer than the direct sum decom- 
position of the underlying Hilbert space T-L into IrlS's of 
<f>. This implication echoes with the findings by Blume- 
Kohout et al.fiJ/ that the essential geometric structure 
underlying all noiseless subsystems, decoherence-free sub- 
spaces, pointer bases and quantum error- correcting codes 
on finite- dimensional quantum systems is an isometry to 
fixed points of certain unital quantum operations. 

IV. DISCUSSIONS 

In summary, using elementary analysis and graph the- 
oretic methods, we completely characterize the structure 
of a unital quantum operation <j> on B(%) provided that 
the Hilbert space T-L is finite-dimensional. In particular, 
the basic building blocks for this kind of unital quan- 
tum operations are those without proper MS induced by 
<j> in the sense that no convex cone formed by the set 
of all positive semidefinite operators acting on a proper 
subspace of H is fixed by </>. We further show that al- 
though the direct sum decomposition of Ti into IrlS's of 
<f> need not be unique, the number of IrlS's and the dimen- 
sion of each IrlS's are unique up to permutation of these 
IrlS's. Using this structural characterization, we solve 
three interesting quantum information processing prob- 
lems, namely, to show a reduction theorem for various 
information-theoretic capacities of a finite-dimensional 
unital quantum channel, to find all the fixed states of 
<f>, and to give a necessary and sufficient condition for 
the preservation of measurement statistics of a POVM 
measurement by <p. 

Interestingly, we find that the problem of completely 
characterizing 4> is reduced to the simultaneous block di- 



agonalization by unitary conjugation of its Kraus oper- 
ators. Note that in actual practice, this simultaneous 
block diagonalization can be done relatively painlessly. 
One possibility is to use a recent computationally stable 
algorithm reported by Maehara and Murota in Ref. [2(| , 
which builds on their earlier works on finding simultane- 
ous invariant subspaces in Refs. [Il|,[22j]. Another pos- 
sibility is to adapt the algorithm of Blume-Kohout et 
al. in Refs. fl2l . Il3j on finding the so-called information- 
preserving structures of a quantum channel. Thinking 
along this direction, the following result of Mendl and 
Wolf [IH may be of interest. They showed that a unital 
quantum operation can always be expressed as the sum 
J2i a iUi ■ U\ where a; el obeying J^i a i = 1- .15] Note 
that this is not an operator-sum representation for can 
be negative. It is instructive to see how to efficiently con- 
vert an operation-sum representation of a unital quantum 
operation <f> to the above affine sum of unitary conjuga- 
tions in a computationally stable manner. It is also in- 
structive to see if Kribs' [l| and Arias et a/.'s [2j necessary 
and sufficient condition on fixed points of a unital quan- 
tum operation restated as Theorem [1] here can be mod- 
ified to cover the case of the affine sum representation 
of Mendl and Wolf. Since finding simultaneous diagonal 
blocks of a set of unitary matrices is much simpler prob- 
lem, solving the above two problems may improve our 
computational efficiency and accuracy in finding IrlS's of 

4>. 
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